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Abstract 

Liouville's theorem, based on the Hamiltonian flow (micro-canonical ensemble) for a many parti- 
cle system, indicates that the (stationary) equilibrium probability distribution is a function of the 
Hamiltonian. A canonical ensemble corresponds to a micro-canonical one at thermodynamic limit. 
On the contrary the dynamics of a single Brownian particle (BP) being explicitly non-Hamiltonian 
with a force and damping term in it and at the other extreme to thermodynamic limit admits the 
Maxwell-distribution (MD) for its velocity and Boltmann-distribution (BD) for positions (when 
in a potential). This is due to the fluctuation-dissipation relation (FDR), as was first introduced 
by Einstein, which forces the Maxwell distribution to the Brownian particles. For a structureless 
BP, that, this theory works is an experimentally verified fact over a century now. Considering a 
structured Brownian particle we will show that the BD and MD fails to ensure equilibrium. We will 
derive a generalized FDR on the basis of the demand of zero current on inhomogeneous space. Our 
FDR and resulting generalized equilibrium distributions recover the standard ones at appropriate 
limits. 
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A small system (mesoscopic system) like a Brownian particle (BP), which is far away 
from thermodynamic limit, can be studied within the framework of Langevin equation 
Langevin equation is Newton's equation of motion for such systems with a dissipative term 
and a force term in it which represent the average interactions of the system with the heat- 
bath. For the system to equilibrate, its energy will, on average (ideally over some span 
of space and time), become a constant and that requires a balance between the energy 
gain from the bath and that dissipated away to the bath. The FDR ensures this energy 
balance. As was first shown by Einstein {2] in some equivalent form, the strength of the 
stochastic force a/2TXbT on a BP (of considerably bigger size than the molecules of the heat 
bath its immersed in) is related to its damping constant T. Such a particle will equilibrate 
with the heat-bath at temperature T with a Maxwellian velocity distribution e - mv2 / 2k BT 



The reason being, the equipartition of energy [4| is inbuilt in the FDR and equipartition 
is a consequence of the Maxwell-distribution of velocity. In the above expressions ks is 
Boltzmann constant, v is the velocity and m is the mass of the BP. Note that, the Brownian 
particle being an object with no structural fluctuations and the heat-bath being homogeneous 
everywhere, the T is a constant. 

Fluctuation dissipation relation in general is an expression of t 
of the system in terms of the correlations of the stochastic forces 
which was given for a conventional Brownian particle considers a constant damping constant 
for the particle where the underlying reality being the BP is large compared to the size of 
the bath degrees of freedom. For Brownian particles of the same dimensions as that of the 
bath degrees of freedom, a classic example has been given by Kubo [8y considering a time 
dependent dissipation for a BP of comparable size of the bath degrees of freedom. The FDR 
is not only useful as an equilibrium condition, its the most widely used instrument to look 
into non-equilirium states and has been proposed in a number of ways and forms |s-12]. In 
the present letter, with a simple model system characterized by a space dependent damping 
r(a;), we will show that a standard form of FDR and the resulting MB distribution does 
not ensure the most important ingredient of equilibrium i.e. zero current (consequently zero 
entropy production) over the associated inhomogeneous phase space of the system. We will 
work out a generalization of the FDR to get a distribution that ensures required zero current 
in equilibrium. This will immediately mean a generalization of the Maxwell-Boltzmann 
distribution which we recover in the standard form over a spatial average. We emphasize on 



re dissipation (response) 
111- Einstein's FDR, 



the fact that in what follows we will be dealing with a system with absolutely no external 
force working on it and the system is immersed in a heat-bath which is uniform everywhere 
i.e. makes the system see the same temperature everywhere in it when in equilibrium. 

The results that we present in this paper are due to coordinate dependence of the dis- 
sipation of small (mesoscopic) structured objects. We will model such a system with no 
external or global field (force) acting on it. Not to tamper with the spatial homogeneity 
of the heat-bath, we will consider a conformation or internal coordinate dependence of the 
damping constant in our model. This is the most important step of the modelling, because, 
an inhomogeneous bath is inherently a non-equilibrium situation. We also have another 
goal in modeling the system in this particular way. As there cannot be any current over an 
inhomogeneous space in equilibrium, there is no reason why there should not be a current 
over a homogeneous space in equilibrium. Any uniform motion in a homogeneous space does 
not produce or destroy anything (including entropy). So, equilibrium will not be disturbed. 
This point is particularly miss-understood by mixing up non-equilibrium steady states with 
a current in the presence of a gradient with the possibility of a current over homogeneous 
space under equilibrium. A recent work of present author shows a particular example of 



such an equilibrium current in a similar model 



13]. We will show that the condition of no 



current in the conformation space practically means the possibility of a uniform current in 



13| 



the bath space and this would be a general proof of the special case mentioned in 
Our model is of the following form 

ri(z)ii = + y^D^) Vl (t) (1) 

T 2 (z)x 2 = + y/2D&jri2(t). (2) 

OX 2 

We consider two particles at the positions x± and x 2 interacting through some potential U (z). 
Moreover, we consider that the damping constant of the particles ^i(z) and T 2 (z) depend 
on the conformation (or internal coordinate) of the system z = X\ — x 2 . The Gaussian 
white noise, that the individual particles experience, have their strengths dependent on the 
conformation of the system because the damping constants are a function of z and that is a 
lesson learnt from the Einstein's FDR. As usual, for a Gaussian white noises, < f]i(t) >= 
and < 771(^1)772 (^2) >= 5i 2 5(ti — t 2 ). We will find out an explicit form for these -D,s from 
the demand of the equilibrium of our system with the heat-bath. Note that, this kind of 
a stochastic noise is actually different from the generally considered multiplicative noise in 



the sense that, the present one does not make the bath space inhomogeneous jl^ . 

Before going into the analysis of our model, let us have a look at the meaning of this con- 
formation dependent damping considered. For the sake of simplicity, consider the simplest 
thing that Stoke's law tells us the damping constant for a spherical body of radius R to be 
proportional to radius R. Thus, R being the internal coordinate defining the conformation 
of the sphere, there is clearly a conformation dependence of the damping constant. Two 
such spheres of radii R\ and R2 being immersed in the same heat-bath will see the same 
temperature (in equilibrium) but different damping constant and different stochastic force 
strengths. In fact, if a sphere changes its radius slowly enough from Ri to R2 being driven 
by the equilibrium fluctuations of the bath, its damping constant and the stochastic force 
strength will also change with the radius while temperature remains the same. The slowness 
condition is required to make sure that equilibrium prevails as the system moves around on 
its conformation space. In our model we have generalized such a system. 

Considering the centre of mass (CM) of the system x = (xi + x 2 )/2, we can recast the 
above equations in the form 

r 1 (z) + r 2 (z)du(z) 



TAz^z) dz 



+ Ut) (3) 



X ~ 2T l (z)T 2 (z) dz 1 j 

where ^ = 772 (*) v ra(z) and = Vi(t) v Ti{z) + ^(t) v r2(z) ■ The cor- 
relations of the £'s can be easily computed using those of 77's to give us < ^2(^1)^(^2) >= 
T^Sfa - h) + g^tx - t 2 ) and < Uh)Ut2) >= - t 2 ) + jggj^ti - h), 

where < £ z (t) >=< £ x (t) >= 0. There also exist cross correlations between £ 2 and ^ 2 which 
we do not require in the present context because of complete decoupling of the conformation 
dynamics from the CM one. 

Let us concentrate on Eq.3, which is equivalent to the motion of a Brownian particle 
in the z-space where the damping depends on the coordinate. The average velocity of the 
particle in the inhomogeneous (due to the presence of U(z) and T(z)s) conformation space 
is given by 

where P{z) is the equilibrium probability distribution of the system which is independent of 
time and the integration limits are at two extreme ends in the z-space where P(z) vanishes. 
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Note the most important point here, that, the standard equilibrium probability distribution 
of Boltzmann form i.e. P(z) = Ne~ u ^ kBT (N is a normalization constant) would surely 
not make the average current vanish for arbitrary F(z)s. The average current identically 
vanishes with Boltzmann-distribution for constant Ts which is a standard result. 

Taking insights from the constant T case, let us impose a sufficient condition for the 
average current to vanish i.e. 

dP(z) 



dz 



CP(z)((z)F(z), (6) 



where C(z) = r ^y^S 2 ,\ which has the dimension of inverse I\ C is a dimension fixing 

^ v i ri(z)r 2 (z) ' ° 

constant and F(z) is the force — dU ^ . This sufficient condition immediately results in a 
probability distribution function 

P(z)=Nexp(C J dz((z)F(z)). (7) 

This distribution function should take the form of Boltzmann distribution at constant T 
corresponding to which ((z) = (const, if the dimension fixing constant is of the form C = 
(constksT and we would keep this in mind. The Fokker-Planck (FP) equation corresponding 
to the Langevin dynamics (Eq.3) can be readily written down as 



dP(z,t) d_ 
dt " ~dz 



C{z)F{z)P{z,t) - 



2£> 1 (z) 2D 2 (z) 



r,(zr r 2 ( z y 



§- z P (^))- (8) 



Now, if we are to produce a stationary probability distribution of the form that in general 
makes < z > vanish, the square bracket term in the diffusion part of the FP equation 
above has to be independent of conformations which is easily achieved if D^z) oc T^z) 2 . 
The proportionality constant is easy to find out keeping in mind the condition already 
mentioned, that, the distribution has to converge to a Boltzmann one for a constant T. The 
proportionality constant comes out to be ksT / < T(z) >. Thus, the generalized probability 
distribution can be written down in the form 

„/ x »r f dz((z)F(z) 

P M = j Vexp ^i±±l , (9) 

where < ((z) >= ^[fz^K^iz)^ • ^ e no ^ e on P assm S t na ^ ^ ne general condition of a 
vanishing < z > implies that the < x > is not in general zero and this is the result showing 
that a current in homogeneous bath space is possible at equilibrium. 



We would generalize our analysis at this point, keeping in mind that we are effectively 
dealing with a single particle scenario in the form of Eq.3 and having known the generalized 
structure of the FDR on the basis of our model, let us apply it to a general single particle 
case which is not over-damped. To that end, we consider the general Langevin dynamics of 
a particle as 

mv(x) = -mT(x)v(x) - F(x) + T(x)J < ^ B > V (t) (10) 

where x denotes the relevant internal coordinate and m is the corresponding reduced mass. 
The above mentioned equation can be obtained from our dimer model easily by setting 
Ti(z) = T 2 (z) to get rid of the involvement of the centre of mass velocity term in this 
equation and we adopt this method for the sake of simplicity. In such a situation the centre 
of mass would undergo a conventional Brownian motion. One can immediately find out 
the expression for the distribution function from the Generalized FP equation which would 
obviously come out to be 

P(v, x) = iVexp ( 2 ^J )k J ) exp ( dx> > ), (11) 

and the analogy is clear with our actual two-body overdamped system. Important point to 
note is, the distribution function does not factorize into position and velocity parts and that 
is expected of a noise having spatial dependence. The normalization constant is given by 
AT" 1 = / dx exp ( f ^ dx' ) / dv(x) exp (-^if^U^)- The global average kinetic 

energy is going to be 

mJV f , . f x < V(x) > F(x'), f , . . . x2 . m< T(x) > v(xf . 

where P(x) is obtained by integrating out v from P(v, x). In the over-damped limit, we can 
easily work out the above mentioned Boltzmann distribution for the position which would 
make <f>=<i>=0 following the line of our previous analysis of the dimer. This is 
an obvious result beyond over-damped regime as well, because, the velocity distribution has 
explicit zero mean everywhere in space whereas the distribution width is a function of space. 

One can also work out a local time average of the energy given to the system by the bath, 
considering the integral of the form 



E 

-'-'in 



^ dt < £(x,t)v(x,t) >, (13) 



where £(x,t) is the stochastic force and <> indicates here a noise-average. From the direct 
solution of the Eq.10, v(x,t) can be found out and v(x,t)/2 is the average velocity within 
the infinitesimal interval dt. Considering the average velocity in this form does away with 
any memory in the velocity and only the velocity developed from zero initial value is taken 
into consideration. The factor of half comes from the approximation of an effectively linear 
dependence of the velocity over time dt under the action of a constant force within this 
interval. This average local supply of energy to the system more explicitly is going to be 

= mY{x)v{xf = E diss , (14) 

where the right hand side is the average energy lost at x with the local average of squared 
velocity v (x) 2 . This energy balance is crucial and ensures equilibrium over the mesoscopic 
time r (bigger than the local correlation time of the noise). This relation implies a local 
conservation of kinetic energy which can be checked easily by multiplying Eq.10 by v(x) and 
taking average over the local velocity. The first and last terms on the right hand side would 
cancel together and the middle term will vanish because < v(x) >= 0, which will mean the 
kinetic energy is locally conserved. 

To conclude, let us have a look at the most important implications of this generalization 
of the FDR and the following distribution. The probability distribution that we obtained is 
the essential finer structure within the MB distribution (that holds globally) which has to 
make the average current vanish over inhomogeneous space in equilibrium. Of course, for 
the average current to vanish, the system has to see its whole conformation space and the 
MB distribution is recovered over that time and that is why we recognize our distribution 
as the internal structure of MB distribution. The importance of our present result is in 
revealing that the barrier overcoming processes in equilibrium, in general, is not governed 
by the potential energy alone as is implied by the standard Boltzmann distribution. Rather, 
the damping actually plays a role in this barrier overcoming. This may prove to be a 
very important step in understanding barrier overcoming processes in polymers, proteins 
(resolving Levinthal's paradox) where the conformation dependence of damping is a reality 
in all likelihood. 

The other important result of our present analysis is the general proof of a uniform 
transport through the heat-bath under equilibrium conditions when the symmetry of the 
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system is broken i.e. Ti ^ T2- This demonstrates Newton's first law in the context of 
equilibrium with a heat-bath that in general causes damping. Its important to note that, 
the present analysis is potentially opening up of a new direction in the understanding of 
complex mesoscopic systems and can be extended to a very large extent. An experimental 
verification of the finer structure of the probability distribution can possibly be probed by 
using polymers or proteins and we hope for such results to come up in future. Considering 
that the expression for the potential U(z) is invertible one can use the potential as the 
variable to write the probability distribution for position as P(z) = Ne~®( u ^ kBT where 
the C[Z(U)]/ < C[Z(U)] >= d$(U)/dU and the force is as usual F = -dU/dZ. This 
expression has the same structure of the Boltzmann distribution. Since, at present we have 
a good idea about what possibly are the forces involved and their corresponding potentials, 
from a distribution that deviates from standard BD, one can guess 3>(U) from curve fitting 
and that would give an idea as to what is the form of the conformation dependent damping 
involved in these cases. 

Another important point to be mentioned is, Eq.10 which is presented here for a homo- 
geneous bath by considering internal coordinate and reduced mass appears to generally hold 
for any coordinate and individual mass of particles. But, this amounts to considering an 
inhomogeneous space in the sense that the damping constant and the external field present 
are making the space inhomogeneous. Through, the form of the present FDR, however, the 
temperature offered by the bath remains the same everywhere. This is a fictitious situation 
where damping is determined by other factors than the bath alone and there exists an FDR 
relating strength of the stochastic force to the damping. So, no FDR of the present form 
should generally be extended into such an inhomogeneous space to account for equilibrium. 



I. DISCUSSION ON ITO VS STRATONOVICH DILEMMA IN RELATION TO 
OUR GENERALIZATION 



We will start this discussion by referring to a recent paper by Lau and Lubensky [15] which 
gives a beautifully concise account of these conventions for tackling stochastic differential 
equations. At the onset we will like to mention that, in our generalization of the FDR 
this controversy over the choice of Ito vs Stratonovich does become irrelevant because of 
the fact that the stochastic force actually appearing in the integrals is no longer position 
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dependent. Note that, the position dependent stochastic force does not only give rise to 
the problem of Ito vs Stratonovich, it also causes a common problem of drift current in 
both the conventions Ito and Stratonovich which is because of the choice of the stochastic 
force to ensure the canonical distribution as an equilibrium distribution of the system at 
large times. This trouble comes alike in Ito convention and Stratonovish convention arising 
from the second moment of stochastic force. This ambiguity is also overcome through our 
method. Let us also mention in the onset that in what follows we will always consider 
that the spatial dependence of the stochastic force is on the internal coordinate. Otherwise, 
a spatial dependence of the stochastic force makes the bath inhomogeneous and that is 
associated with many other problems like sustenance of this inhomogeneity in the presence 
of bath related drift and the entropy production associated with those processes etc which 
are not generally addressed while claiming equilibrium on such systems, but, are genuine 
issues related to equilibrium. 

As has been shown in the ref. 15) (and in numerous other generalization of Stratonovich 



convention), a stochastic differential equation of the form 

% = -n^ + g(xHt), (15) 

where t](t) is a white noise. This form of an equation is considered not enough for the 
system to go to equilibrium at large times, one has to add one more drift term in a generalized 
Stratonovich approach. The term one adds (following ref. [2]) is f% = (l—a)g(x)g'(x) (where 
prime indicates derivative with respect to x). When a = the convention is Ito, a = 1/2 
its Stratonovich and a any other value gives us a generalized Stratonovich convention. Its 
interesting to note that, even in Ito (a = 0), there is a need for an additional drift term 
which is g(x)g'(x) and the sole origin of this term is the demand of a canonical distribution 
in equilibrium for the system to remain thermodynamically correct at the extreme other end 
of thermodynamic limit! This drift current is a result of considering the stochastic force as 
g 2 (x) = 2T(x)kBT to ensure canonical distribution and this problem is not there if g(x) is 
a constant resulting in an effectively constant-width stochastic contribution to the position 
fluctuations. 

There are two crucial ingredients required for ensuring equilibrium at large times in the 
above mentioned methods, 1. the current is made zero over inhomogeneous space to avoid 
entropy production and maintaining microscopic detailed balance and 2. the probability 



distribution has to have a canonical form for equilibrium. The price paid is in requiring 
bath driven drift which is done away with by adding an extra drift term in the model and 
thereby essentially requiring the noise to be in a form consisting of a deterministic and 
a stochastic term. The second point is even more interesting, it shows that the relation 
required for equilibrium i.e. g 2 (x) = 2T{x)ksT is local and the system with inhomogeneity 
over space is not required to see the whole space over which it will actually equilibrate to get 
the strength of the stochastic force fixed for equilibrium. We argue that, a simple Brownian 
particle when confined by a potential gets a canonical equilibrium distribution because of 
the fact that it sees a uniform bath and thus the locality is not an issue. But, for a particle 
which sees a non-uniform heat bath and has to equilibrate over a finite region being bound 
to a potential, must see the whole space available to it to equilibrate. 

In our modification of the expression of the stochastic force for the equilibrium of the 
system we have relaxed the above mentioned condition of canonical probability distribution 
keeping the other one of requiring zero current there. We have done so based on the convic- 
tion that the thermodynamics is not exactly the same away from thermodynamic limit. The 
most crucial advantage of this method is that the stochastic contribution to the position 
fluctuations (at the over-damped limit) of the system is a constant and has the strength 
2k B T/ < F(x) > where < T(x) >= j dxT(x)P(x) where P(x) is the probability of x. This 
is a strength of the force which depends on the global parameter <F(x) > strictly requiring 
the system to see the whole inhomogeneous space available to fix the equilibrium strength 
of the stochastic force. In this way, there is no need of a deterministic part added to the 
stochastic force to get rid of spurious drifts and there is no controversy of Ito vs Stratonovich 
and not even Ito and Stratonovich. For the model with inertial terms present, it can also 
be converted to one with an effective stochastic force with constant strength by a change 
of variable from v to v/F(x). The Probability distribution is in a form which resembles 
Boltzmann form more closely because of not having any space dependent amplitude to the 
exponential part of the probability distribution. This constant amplitude form of the prob- 
ability distribution is required for the symmetry j dxP(x)f Lg = J dxP(x)gLf to exist, 
where / and g are any doubly differentiable functions and L = —T(x)^--j^ + g{x)-^. This 
symmetry ensures detailed balance of the system and to maintain this the spurious drift 
terms are actually added in Ito and Stratonovich approach. The probability distribution of 
ours will converge to canonical distribution when internal degrees of freedom, over which the 
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inhomogeneity exists, are integrated out. Important to note that the fixing of the constant 
in the strength of the equilibrium stochastic force to < T(x) > is required by the recov- 
ery of equipartition of energy over homogeneous space as is shown in Eq.12. Now, which 
method represents what actually happens in nature has to be determined by experiments 
given the apparent scarcity of information for writing down a dynamics for the evolution of 
the stochastic force to its equilibrium value. 
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